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Abstract 

The dynamic mechanical response of mesoporous sihca with coated inner surfaces confining the 
glass forming liquid salol is measured as a function of temperature and frequency (1-100 Hz) for 
various pore sizes (2.4-7.5 nm). Compared to former results on natural pores, a distinct acceleration 
of dynamics due to the removal of surface-related retardation of molecular dynamics is found now, 
which can be fitted by a homogeneous relaxation using an unmodified Vogel-Fulcher-Tammann 
relation. This lubrication effect leads to a stronger decrease of the glass transition temperature Tg 
with decreasing pore size. The present data allow to quantify and separate competing side effects as 
surface bondings and negative pressure from the pure confinement induced acceleration of molecular 
dynamics with decreasing pore size. We analyze the dynamic elastic susceptibility data in terms 
of a recently proposed procedure [C. Dalle- Ferrier et al., Phys. Rev. E 76, 041510 (2007)], which 
relates the number Ncorr of molecules that are dynamically correlated to the three-point dynamic 
susceptibility. Assuming that Ncorr is representative of the size ^ of a dynamically correlated region, 
i.e. Ncorr ^ we find that ^ grows when approaching the glass transition to about 3 nm, a size 
which was also reported previously using other techniques. 
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I. INTRODUCTION 



Glass- forming materials have been produced by mankind for more than 6000 years. 
Despite several decades of intense research the transition of a liquid into its glassy state is 
still lacking a universal theory explaining both the increase of viscosity r] and molecular 
relaxation rates by 14 orders of magnitude |l, 2| without creating any long range order. A 
widely used explanation which goes back to Adam and Gibbs |I3] is based on the assumption 
of cooperative rearrangement of molecules ("cooperative rearranging regions, CRR"), 
forming compact clusters of a typical size ^. Such a subsystem of molecules can rearrange 
into another configuration independently of its environment. 

The size of these groups of molecules is considered to grow to some nm as Tg is approached 
{4]. E.g. random first order transition theory of glasses predicts ^ = ro0.51 {In—)^. 

At Tg where r 100 s one obtains for typical values of tq ~ 10~^^ s and Tq ~ 1 nm, 
iiTg) ^ 5 nm. Measuring such dynamic heterogeneities is one of the most important but 
at the same time difficult issues in the field of glass formation. Several experimental setups 
hke specific heat spectroscopy_Q, Q|, multidimensional NMR js-ioi and non-hnear dynamic 
susceptibility measurements [11], etc. have been used to determine a possible growing 
length scale accompanying the glass transition. All these results agree in the fact, that the 
obtained cooperative regions are of the order of several nm near Tg, and are displaying 
a weak temperature dependence [l2|- Recent computer simulations also confirm a rather 
modest growth of the dynamically correlated regions when approaching th e g lass transition 



from above [13|. Very recently in a groundbreaking work Biroli et al. found direct 
evidence for a growing dynamical length in supercooled liquids by applying inhomogeneous 
mode-coupling theory. Based on this theory these authors have recently developed a 
method to quantify the size of the correlated regions [ll| by analyzing three-point dynamic 
susceptibilites. For a large number of supercooled liquids they could indeed confirm growing 
of the correlated regions when approaching the glass transition (isl . 

However, direct measurements of ^ are extremely difficult, sometimes even impossible 
and one has to resort to indirect investigations. A widely used approach is confining the 
glass-forming liquid spatially, either in the form of thin films or by using mesoporous 
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host matrices for confinement. If the size d of the confinement is finite the cooperatively 
rearranging regions cannot grow beyond any bound, becoming saturated at = (i < oo. 
This should lead to a confinement induced acceleration of the dynamics resulting in a 
downshift of the glass transition, and even impeding it |l6| as > c/. 

Since the pioneer work of Jackson and McKenna |r?| in 1991, uncovering a reduction 
of the glass transition temperature AT^ cx 1/d in confinement of size (i, a variety of 
confinement geometries and experimental methods have been used. Both weak and strong 
glass forming liquids, showing strong and weak interaction with the 2D or 3D confinement 
media were investigated. Single, double, even multiple transitions have been observed. 
Extensive topical overviews are found in Refs. Q, Sandy. 

The abundance and diversity of experimental findings shows that an accurate discussion 
of side effects in discussing results of glass forming liquids in confinement is essential. 



Negative pressure due to mismatching thermal expansion coefficients o 



liquid and confining 



matrix is such a side effect. It was discussed by various authors j21,j22| and sometimes 
even made responsible for the whole downshift of Tg in confinement [21'|. Being true, it 
would disprove the idea of a growing length scale of cooperativity. In a former paper (23| 
the authors have determined negative pressure effects for salol in natural un-coated pores of 
size 7.5 to 2.6 nm from high resolution thermal expansion measurements. An upper bound 
for the contribution of negative pressure to the total downshift of Tg of ~ 30% was found. 

A second, and much larger effect on the glass transition of liquids in confinement arises 
from the interaction of the molecules with the huge inner surface of confining host matrices 
which can take values up to 600 w? /g (see Tab. Ij). Confined liquids tend to form H-bonds 
with the hydrophilic pore surface, which leads to an immobile surface layer of molecules 
and a retarded relaxation behavior at the glass transition. 

In recent dynamic elastic measurements of salol filled into matrices of Vycor and 

Gelsil with natural untreated pores we studied this competition between surface induced 
slowing down and confinement induced acceleration of the dynamics. Here we present new 
results of salol confined in mesoporous Vycor and Gelsil with silanated pore surfaces. The 
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results clearly show that silanation reraoves the liquid-surface interaction, leading to an 
enhancement of the molecular dynamics in the pores, resulting in a stronger downshift 
of Tg as compared to the uncoated pores. These findings allow to separate the surface 
effect from confinement induced acceleration, and a simultaneous quantitative statement 
about negative pressure within one and the same measurement technique and confinement 
geometry. A comparison with previous dielectric spectroscopy data j25-27| of salol in 
mesoporous matrices yields excellent agreement. 



Using the recently proposed method of Berthier et al. 



we determine the number 



NcOrr of dynamically correlated molecules as a function of temperature and pore size. 



II. EXPERIMENTAL 



A. Sample preparation 



Vycor by Corning Inc., NY is produced via phase separation within a Na20-B203-Si02 



melt, and subsequent acid leaching [28|, which leaves a 98% pure Si02 skeleton containing 
interconnected pores of random length, direction and density. A very narrow pore size 
distribution and an average ratio of pore length / over pore diameter d of l/d ^^4.35 is 
reported ^ . Gelsil samples result from a sol-gel process [30] and consist of randomly 
formed pure fused Si02 monodisperse spheres (3l|, touching and penetrating each other, 
resulting in a mesoporous structure with a rather broad distribution of pore diameters, as 
N2-sorption results showed. Samples are cut and sanded to gain the required orthogonal 
shapes, and cleaned in a 30% H2O2 solution at 363 K for 24h. Drying is done under high 
vacuum (10~^ mbar) at 393 K for another 24h. 



In order to deactivate inner surfaces, 0H~ groups are replaced by OSi(CH3)3 trimethylsi- 
lyl groups via exposing cleaned samples to gaseous hexamethyldisilazane (HDMS, from 
Sigma Aldrich, purity 99.9%) in a closed vessel at 330 K for 24h. Afterwards samples are 
again evacuated for 24h at 300 K. In order to check pore geometry and pore radius after 
silanation, mesoporous samples have been tested via N2-sorption and BET/BJH analysis 
of the individual desorption isotherms {32! experiments before and after HDMS treatment. 
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Results are shown in Tab. ID Whereas adsorption data on Vycor does not show a significant 
change of pore geometry, in Gelsil samples surface area and pore volume were found to de- 
crease strongly due to silanation. The effect of silanation on the pore diameter is ^ 0.2 nm, 
corresponding to a HDMS layer thickness of ~ 0.1 nm. For comparison Kremer et al. [25| 
estimated the thickness of the silan layer from analysis of dielectric strength data as 0.38 nm. 

Filling with the fragile low molecular weight glass former jsSi phenylsalicylate [salol, 
C13H10O3, T^(bulk)= 316 K] is done by capillary wetting at 333 K for 12h. Filling fractions 
/ are evaluated via weighing clean and filled samples. 

Some characterizing parameters for the glass-freezing behavior of bulk salol are: fragility 
index m - which usually varies from m=17 for strong glass formers to m=150 for fragile 
ones - takes for salol [34] the value m=73, T^(bulk)= 220 K (glass transition temperature 



defined at r=100 s) and Typ=175 K (35| (Vogel-Fulcher temperature). The volume of a 
salol molecule is estimated (36| as 0.282 nm^ corresponding to a mean diameter of about 
0.8 nm. 



B. Dynamic mechanical analysis 



The dynamic mechanical response of samples of a typical size of 1 x 2 x 8 mm^ in three 
point bending (3PB) mode yields the real and imaginary part of the complex Young's mo- 



dulus Y* = -{-iV within a frequency range of 0.01 to 
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00 Hz applying static and dynamic 



371 and 



38l . Temperature is controlled 



forces up to 9.8 N. For further details see Refs. |23|, 
by gaseous nitrogen flow from 120 K to RT. DMA analyzers are decoupled from building 
vibrations, and electronics are shielded from a possible interference with the 50 Hz mains 
voltage frequency. The analyzers used are a series 7 DMA and a Diamond DMA, both built 
by Per kin Elmer Inc. 



C. Resonant ultrasound spectroscopy 



Due to contact errors a DMA experiment does not yield absolute values for elastic moduli. 
Therefore resonant ultrasound spectroscopy (RUS) was used to determine bulk and shear 
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TABLE I: N2 adsorption characteristics and elastic moduli of untreated and silanated porous silica 
samples. 





Gelsil 2.6 


Gelsil 5 


Vycor 


untreated 


av. pore diameter (nm) 


2.6 


5.0 


7.5 


surface area (m^/g) 


590 


510 


70 


pore volume (cm^/g) 


0.38 


0.68 


0.21 


porosity (f) 


0.36 


0.54 


0.40 


bulk mod. K (GPa) 


9.6 


3.9 


8.1 


shear mod. G (GPa) 


7.7 


3.3 


6.7 


Young 's mod. Y (GPa) 


18.2 


7.7 


15.8 


silanated 


av. pore diameter (nm) 


2.4 


4.8 


7.5 


surface area (m^/g) 


260 


325 


65 


pore volume (cm^/g) 


0.15 


0.4 


0.19 


porosity (j) 


0.30 


0.49 


0.33 


bulk mod. K (GPa) 


9.6 


3.3 


9.1 


shear mod. G (GPa) 


9.0 


3.9 


9.0 


Young's mod. Y (GPa) 


20.6 


8.9 


20.3 



moduli of both natural and silanated mesoporous samples at room temperature. Orthogonal 
parallelepipeds of about 2.9 x 3.0 x 3.1 mm^ were used to gather resonance spectra from 
50 to 1100 kHz. For each sample 25 peaks resulting from excited resonant eigenmodes and 
corresponding overtones then were fitted via a Lagrangian minimization routine gaining 
bulk modulus K and shear modulus G (see Tab. Ij) with an accuracy of less than 1%. For 
further experimental details see e. g. Ref. 39. Young's modulus Y was calculated from 
Y = 9KG/ {3K -\-G) and used to calibrate DMA raw data at room temperature (see Figs. [H 
ElandlSj). 
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III. RESULTS AND DISCUSSION 



As an example, the dynamic elastic response (IHz - 100 Hz) of salol in Gelsil with 
silanated pores of 4.8 nm diameter is shown in Fig. [TJ The data for 2.4 and 7.5 nm pores 
look very similar and are hence not presented here. 

Fig. |2] gives a comparison between the recent results of salol in natural untreated pores 
{23! and the new data on silanated samples. The most striking feature of the silanated 
samples is the absence of a double peak structure and the shoulder in Y'\T) and Y\T), 
respectively. In untreated samples a significant part of the confined liquid sticks to the pore 
surface due to hydrogen bonding, thus being retarded in relaxation dynamics. This leads to 
a local dependance of relaxation times across the pore section (see Fig. Hj) causing an extra 
glass transition at higher temperature relative to the one of the core molecules. Assuming a 
spatial distribution of Vogel-Fulcher temperatures [40] To(r) = Too -\- k/{R — r -\- Vp), which 
translates via r(r) = tq • exp[E/(T — To(r))] to a distribution of relaxation times we were 
able to fit the stepwise change in Y'(T) and the double-peak-structure in Y"{T). A detailed 



analysis can be found in Ref. 



23 



For silanated surfaces we now find only bulk-like relaxation, i.e. just one step in Y'{T) 
and one single narrow peak in Y"{T). Peaks do shift with measurement frequency, as Fig. [1] 
shows. Since the liquid- surface interaction is removed now, we do not take into account 
any radial dependence of the relaxation time (see Fig. Hj), but use a single homogeneous 
Vogel-Fulcher- Tammann equation 

r(T)=ro.e.p[,^^] (1) 

with the pre-exponential factor ro, the activation energy E and the Vogel-Fulcher tempera- 
ture To^d), depending on the pore diameter d. 



Similar as in our previous work |23|, |2^ a Cole-Davidson relaxation is used to model 
dynamic mechanic susceptibility in terms of (now radially monodisperse) relaxors 

y\^) ^ ^ ^ (2) 
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Gelsil 4.8 nm silanated + salol 
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200 



220 
T(K) 



240 



260 




FIG. 1: Real (a) and imaginary (b) parts of the complex Young's modulus of silanated Gelsil 
(4.8 nm) filled with salol (filling fraction / = 0.87) measured in three point bending geometry 
(Diamond DMA). IHz signal are original data, other signals are offset for sake of clarity. 



with uo = 27rz/, v being the measurement frequency, and the broadening parameter 
Separating real and imaginary part of V —Y' + lY" leads to 

cos [7 • arctan{uoT{T))] 



m 7- 



Y'{T) = 1 

y"(T) 



Ay ■ 
Ay • 



sm[7 ■ arctan{LOT{T))] 



(3a) 
(3b) 



[l+Cj2^(T)2]i 

Eqs. (3a, b) and Eq. (P) are used to fit the data given in Figs. [1] and El We point out 
that the pore radius R — d/2 enters as a fitparameter in the pore size dependence of the 
Vogel-Fulcher temperature To{d). Fits and corresponding parameters are shown in Fig. |3] 
and Tab. |IT1 respectively. 



Silanation causes two main effects on the dynamics of the molecules within the pores: 
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FIG. 2: Comparison of DMA data: thin lines are data for untreated inner surface, thick hues 
are new data for silanated pore surfaces. Data are presented in arbitrary units since in a DMA 
experiment contact losses do not allow to quantitatively compare Y' and Y" signals of two different 
samples. 



1.) For untreated surface, we had to consider pore center relaxation times about two orders 

n 

of magnitude higher [23] than obtained from the present analysis of silanated pores. Such 
an enhancement of mobility of molecules due to the absence of surface blocking in silanated 



pores was also observed in previous studies 



271 . l42| . It probably reflects the fact, that 



the surface blocking of molecular mobility in uncoated pores slows down the dynamics of 
molecules also in the center of the pores. 

2.) Present results and used fit parameters (Tab. |Ill) show that the acceleration of 
the dynamics in silanated pores now leads to a much stronger downshift of Tg with 
decreasing pore size as compared with the results of uncoated pores [23^. It reflects the 
pure confinement effect - i.e. free of surface contributions - which leads to a downshift of Tg 
due to limitation of the correlated regions by the pore diameter. 



Unfortunately there is no unique theory which clearly relates these finite size effects to 
par ameters characterizing the glass transition. However in recent computer simulations 



43| of supercooled polymer films confined between two separated walls of distance d it was 



shown that confinement leads to faster dynamics. The authors parametrized the size depen- 
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dence of the relaxation time as r(T, d) oc exp 



E{d) 



and found similar (i-dependencies 



T-ro(d) 

for the mode-coupling critical temperature Tc{d) and the Vogel-Fulcher temperature T^i^d) 



arguing, that Tg{d) should also follow a similar curve (see Fig. 14 of Ref. |43|). As already 
mentioned above we here use the same dependence of the relaxation time (Eq. ^) to fit 
our data, and obtain a very similar size dependence of Tg{d) and T^i^d) (Fig.Ej) as Varnik et al. 



A. Hunt [44j has calculated finite size effects on the glass transition temperature in glass- 
forming liquids analytically using percolation theory, yielding 

T,{d)=T^"'(l-J-j) (4) 

where tq is the typical distance between molecules, which we approximate as the diameter 
of a salol molecule Tq ^ 0.8 nm. Using a Gaussian distribution of energy barriers, the 
author obtained a value oi c ^ 0.27. A fit of our experimentally determined Tg{d) (Table HB 
with Eq. (ij) yields c ^ 0.13 and T^""^^ ^ 216 K (Fig. Ej), which is in rather good agreement 
with the theoretical estimate. In this picture Tg is reduced by confinement because the 
average barrier height for the molecules in pores is smaller than the so called "blocking" 
barrier, which is responsible for the glass freezing in the bulk. We have also plotted the 
values of To(c/) of Table HIl in Fig. [71 The corresponding fine was drawn using the relation 
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FIG. 4: Modeled relaxation time in untreated and silanated pores of Gelsil 5 from Eqn. ([T]) used in 
Ref. 23 and in Eqs. (3a,b) for fits of data in Fig. |3]3 and e herein. 
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FIG. 5: Relaxation time in pore centers calculated from Eq. ([T]) with corresponding parameters 
from Tab. [Hi Horizontal line shows r = 100 s. Gray lines are relaxation times in untreated pores 
from Ref. 23I. 



Tg — To = in{ioo/To) ^ ^^'^ which is obtained from Eq. ([T]) using the fit parameters 
E = 1750 K and tq = 10"^^ s of Table HI 
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FIG. 6: Shift of glass transition temperature against (pore diameter) ^. Open circles display the 



maximum negative pressure contribution (see chapter III. C. of Ref 



same reference, filled triangles show literature values from Ref. 



21 



23l ). boxes are AT^s from the 



Open triangles are the present 



results, and crosses mark corresponding literature data from Ref. 



25| 



As already mentioned above, Dalle- Ferrier et al. [15] have given an expression for the 
number A^corr,4 of molecules that are dynamically correlated over a time interval of the order 
of r as 



Nccyrr,i\T) 



max. 



ACp V " ""^ dT 



(5) 



where ACp in units of the gas constant R is the excess specific heat of the glass-forming liquid 
at constant pressure ist and x{T, lj) is a suitable dynamic correlation function. Very often 
glass-forming materials are studied by dielectric spectroscopy measurements and therefore 
the dynamic susceptibihty is identified with the dielectric susceptibility. To estimate the 
number of dynamically correlated molecules, we apply two different procedures: In the first 
we are using x{TjU)) := ^'^(o)Zy'{^) (0) "to analyze the data directly, i. e. without 

any fit procedure in between. In the second case we rewrite Eq. (JSj) with Eq. (I2j) yielding 



fdlrvT\ 



(6) 



ir o rii t^i — \ ^ r / \ sin[^ arctan{—)] + — cos[^ arctan{ — )] 

where t(T) is obtamed from fits or the data m Fig.|3|and 7(7) = (i^^-2]i+7/2 ~ 

results from the Cole-Davidson dynamic response function and is the analog to the stretched 
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FIG. 7: Shift of glass transition temperature Tg and Vogel-Fulcher temperature Tq against pore 
diameter d. The points are determined from fitting the experimental data with Eqs. (3a, b) and ([T]), 
where Tg[d) is obtained from T{T,d) = 100 s, as shown in Fig. (|5j). The lines are fits using Eq. (j4|) 



and the relation T. 



9 ^0 — ln{100/To) 



50.6 K, as indicated in the text. 



3.5 
3 

2 
1.5 

200 210 220 230 240 250 260 270 

T(K) 

FIG. 8: Size ^ of dynamically correlated regions. Lines are calculated from Eq. ([6j) using the 
procedure described in the text. Symbols are calculated from DMA data using Eq. ([5]). Inset shows 
^ at T = for various pore sizes. 
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TABLE II: Fit parameters used in Eqs. (3a,b) for fits of Fig. E 





Vycor 


GelsilS 


Gelsil2.6 


^(K) 


1750 


1750 


1750 


Too (K) 


161.5 


158.5 


154.5 


^0 (s) 


10-13 


10-13 


10-13 


7 


0.19 


0.25 


0.17 



exponential p of the KWW-response function used e.g. in Eq. (8) of Ref. |45|. Fig. |8] com- 
pares the temperature dependence of the dynamic correlation length ^ calculated from the 
two different methods and applying the hypothesized relation ^ = 2 (3A^corr,4Ka/o//47r)^^^, 
which results in excellent agreement. We used DSC values [33| of ACp{d) reported for salol 
confined in mesoporous silica in a wide range of pore sizes from 2.6 to 26.4 nm. The in- 
set of Fig. |8] shows that the dynamic correlation length at the glass transition temperature 
(,{Tg) for various pore sizes determined from Eq. ([6]) is rather constant with a mean value of 
about 3 nm. This overall behaviour of (,{Tg^ d) is very similar to calorimetrically determined 
characteristic lengths for salol in confined geometries [46 1. 



IV. CONCLUSIONS 



Results of extensive dynamic mechanical measurements of the glass forming liquid 
salol confined in mesoporous silica with silanated pores are presented. It turns out, that 
silanation can completely remove the liquid- surface interactions. As a result the confinement 
induced acceleration of the dynamics can be measured purely, which leads now to a much 
stronger (as compared to uncoated pores) downshift of the glass transition temperature 
Tg with decreasing pore size due to the hindering of cooperativity. Using the results of 
percolation theory [4^ we have calculated the downshift of Tg with decreasing pore size, 
which fits our data very well. In Ref. 44 it is also shown, that finite size effects are expected 
to set in when the pore size 7ro — lOro yielding 4.8 — 8 nm for salol with Tq ^ 0.8 nm. 
This is in very good agreement with our observations. 
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We have also analyzed our dynamic elastic data obtained for the different pore sizes 



in terms of a newly proposed theory 



Uj which relates the size of dynamically correlated 



regions to the temperature derivative of the dynamical two-point correlation function, 
which in our case can be identified with the dynamic susceptibility Y{uj,T). The results 
clearly show an increase of the dynamic correlation length ^ with decreasing tempe- 
rature approaching ^(T^) ~ 3.2 nm. This value is very similar to the one obtained 
previously [4] from DSC and TMDSC data. For smaller pore sizes (,{T,d) at a given 
temperature shifts to smaller values which is concomitant to the systematic decrease of 
the relaxation time r{T^d) and the resulting downshift of Tg{d). However at the glass 
transition temperature the dynamic correlation length is almost independent of the pore size. 



Unfortunately at present there is no unique theory that relates the relevant parameters 
controlling the confinement effects in glass-forming materials to experimental data. The 
reason for this is that the microscopic mechanism behind the glass transition is still not 
completely understood and more theoretical (e.g. of the type presented in Ref. \^ and 
experimental work is required to close the gap of knowledge and understand confinement 
effects in glass-forming liquids. 
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